Random walks and orthogonal polynomials: some 

challenges 

F. Alberto Grunbaum 

Abstract. The study of several naturally arising "nearest neigh- 
bours" random walks benefits from the study of the associated 
orthogonal polynomials and their orthogonality measure. I con- 
sider extensions of this approach to a larger class of random walks. 
This raises a number of open problems. 
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1. Introduction 

Consider a birth and death process, i.e., a discrete time Markov 
chain on the nonnegative integers, with a one step transition probability 
matrix P. There is then a time-honored way of writing down the n-step 
transition probability matrix P n in terms of the orthogonal polynomials 
associated to P and the spectral measure. This goes back to |KMcG| 
and, as they observe, this is nothing but an application of the spectral 
theorem. One can find some precursors of these powerful ideas, see for 
instance [HI ILR| . In as much as this is such a deep and general result, 
it holds in many setups, such as a nearest neighbours random walk on 
the iVth-roots of unity. In general this representation of P n allows one 
to relate properties of the Markov chain, such as recurrence or other 
limiting behaviour, to properties of the orthogonality measure. 

In the few cases when one can get one's hands on the orthogonal- 
ity measure and the polynomials themselves this gives fairly explicit 
answers to various questions. 
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The two main drawbacks to the applicability of this representation 
(to be recalled below) are: 

a) typically one cannot get either the polynomials or the measure 
explicitly. 

b) the method is restricted to "nearest neighbour" transition prob- 
ability chains that give rise to tridiagonal matrices and thus 
to orthogonal polynomials. 

The challenge that we pose in this paper is very simple: try to 
extend the class of Markov processes whose study can benefit from a 
similar association. 

There is an important collection of papers that study in detail the 
cases where the entries in P depends linearly, quadratically or even ra- 
tionally on the index n. We make no attempt to review these results, 
but we just mention that the linear case involves (associated) Laguerre 
and Meixner polynomials, and the other cases involve associated dual 
Hahn polynomials. For a very small sample of important sources deal- 
ing with this connection see |Chl IvDl IILMV| . 

The plan for this paper is as follows: 

1) we first review briefly the approach of S. Karlin and J. Mc- 
Gregor. This is done in section 2. 

2) we consider a few examples of physically important Markov 
chains that happen to feature rather well known families of 
orthogonal polynomials. This is done in sections 3,4,5. 

3) we propose a way of extending this representation to the case 
of certain Markov chains where the one-step transition prob- 
abilty matrix is not necessarily tridiagonal. For concretness 
we restrict ourselves to the case of pentadiagonal matrices or 
more generally block tridiagonal matrices. This is illustrated 
with some examples. This material takes up sections 6,7,8,9 
and 10. 

4) A number of open problems are mentioned along the way. A 
few more are listed in section 11. 

Both parts 1) and 2) above are well known while the proposal in 3) 
appears to be new. 

2. The Karlin-McGregor representation 

If we have 

P ij3 - = Pr{X{n + 1) = j\X(n) = i} 
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for the 1-step transition probability of our Markov chain, and we put 
Pi = P^i+i, qi+i = Pj+i,j, and r\ = P^ we get for the matrix P, in the 
case of a birth and death process, the expression 

fro Po \ 
qi ri pi 
q 2 r 2 p 2 

\ ■■■/ 

We will assume that Pj > 0, q j+1 > and r,- > for j > 0. We also 
assume Pj + rj + q 3 : — 1 for j > 1 and by putting p + r < 1 we allow for 
the state j = to be an absorbing state (with probability 1 — po — r ). 
Some of these conditions can be relaxed. 

If one introduces the polynomials Qj(x) by the conditions <5-i(0) = 
0, Qo(x) = 1 and using the notation 




we insist on the recursion relation 

FQ(x) = xQ(x) 

one proves the existence of a unique measure ijj(dx) supported in [—1, 1] 
such that 

TTj J Qi(x)Qj(x^(dx) = S i:j 
and one gets the Karlin-McGregor representation formula 

(F")* = ^ j x n Q l {x)Q j {x)ij{dx). 

Many general results can be obtained from the representation formula 
given above. Some of these will be given for certain examples in the 
next three sections. 

Here we just remark that the existence of 

lim (P n )^ 

n^oo 

is equivalent to ijj(dx) having no mass at x — — 1. If this is the case 
this limit is positive exactly when ip(dx) has some mass at x — 1. 

If one notices that Q n (x) is nothing but the determinant of the 
(n + 1) x (n + 1) upper-left corner of the matrix xl — P, divided by the 
factor 

PoPi...p„_i, 
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and one defines the polynomials q n (x) by solving the same three-term 
recursion relation satisfied by the polynomials Q n (x), but with the 
indices shifted by one, and the initial conditions qo(x) = 1, qi{x) = 
(x — ri)/pi then it is clear that the (0, 0) entry of the matrix 

(xl -F)- 1 

should be given, except by the constant po, by the limit of the ratio 

q n -i(x)/Q n {x). 

On the other hand the same spectral theorem alluded to above 
establishes an intimate relation between 

h*m n _ 0O g n _i(a;)/Q n (x) 

and 

r 1 #(a) 

We will see in some of the examples a probabilistic interpretation 
for the right hand side of the expression above in terms of generating 
functions. 

The same connection with orthogonal polynomials holds in the case 
of a birth and death process with continuous time, and this has been 
extensively described in the literature. The discrete time situation 
discussed above is enough to illustrate the power of this method. 



3. The Ehrenfest urn model 



Consider the case of a Markov chain in the finite state space 0, 1, 2, 2N 
where the matrix P given by 



/0 1 


2N 



\ 



2N-1 
2N 

-2- 

2N u 



2N-2 
2N 





V 





2N 
2N 



1 

2N 

0/ 



This situation arises in a model introduced by P. and T. Ehrenfest, 
see [E], in an effort to illustrate the issue that irreversibility and re- 
currence can coexist. The background here is, of course, the famous 
if -theorem of L. Boltzmann. 
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For a more detailed discussion of the model see (EJEJ. This model 
has also been considered in dealing with a quantum mechanical version 
of a discrete harmonic oscillator by Schroedinger himself, see [SchKj. 

In this case the corresponding orthogonal polynomials (on a finite 
set) can be given explicitly. Consider the so called Krawtchouk poly- 
nomials, given by means of the (truncated) Gauss series 

a, b \ ^ {a) n {b) n 



3*1 



; z 



E 



n\(c) n 



with 



(a) n = a(a + 1) . . . (a + n — 1) 
The polynomials are given by 

Kit?) = &{~-in> 

x = 0,1,..., 2JV; i = 0,l, 

Observe that 

K (x) = l,K i (2N) = (-l) i . 
The orthogonality measure is read off from 

'2N 

2N 



Wo 



2N 



x=0 



>27v - (-2Ny ij - 1 13 



< i,j < 2N. 



These polynomials satisfy the second order difference equation 

— xKAx) 



\{2N - i)K i+1 {x) - \2NK t {x) + '-K^x) 



2 y ' ' K ' 2 

and this has the consequence that 

/o 

i 

27V 



1 



2 
2N 



2N-1 
2N 





\ 



2N-2 
2N 





^- 

VJ O AT 



\ 



m o j 

2N u / 



2N 



( K (x) \ 
KAx) 



x \ 

n) 



( K (x) \ 
Kt(x) 

\K 2N (x)J 



\K 2N (x)J 



any time that x is one of the values x = 0, 1, . . . , 2N. This means that 
the eigenvalues of the matrix P above are given by the values of 1 
at these values of x, i.e., 



X 

N 



i i 
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and that the corresponding eigenvectors are the values of [K (x), K 1 (x) , . . . , K 2 n(x)] t 
at these values of x. 

Since the matrix P above is the one step transition probability ma- 
trix for our urn model we conclude that 

2N I 

x=0 

We can use these expressions to rederive some results given in [K]. 
We have 

2N 

and the "generating function" for these probabilities, defined by 

oo 
n=0 

becomes 



x \ n \ x 



2 2N 



2N , r 

U ^ = ^ Ml -z) +xz 2 2N ' 

In particular £/(l) = oo and then the familiar "renewal equation", see 
[F] given by 

U(z) = F(z)U(z) + 1 

where F(z) is the generating function for the probabilities f n of return- 
ing from state to state for the first time in n steps 



F(Z) = ^fn 



n=0 

gives 

™ = 1 ~W) 

Therefore we have F(l) = 1, indicating that one returns to state with 
probability one in finite time. 

These results allow us to compute the expected time to return to 
state 0. This expected value is given by F'(l), and we have 

F'(z) = ™ 
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Since 

'2N 



2/V 



U'(z) 



Wa 



N(N-x) \x 



{N{l-z) + xz) 2 2 2N 



x=0 

we get F'(l) = 2 2N . The same method shows that any state i = 
0, . . . , 2N is also recurrent and that the expected time to return to it 
is given by 

2 2N 



2N 



The moral of this story is clear: if i = or 2N, or if i is close to these 
values, i.e., we start from a state where most balls are in one urn it will 
take on average a huge amount of time to get back to this state. 

If on the other hand i = N, i.e., we are starting from a very balanced 
state, then we will (on average) return to this state fairly soon. Thus 
we see how the issues of irreversibility and recurrence are rather subtle. 



In a very precise sense these polynomials are discrete analogs of 
those of Hermite in the case of the real line. For interesting material 
regarding this section the reader should consult [A]. 



4. A Chebyshev type example 

The example below illustrates nicely how certain recurrence prop- 
erties of the process are related to the presence of point masses in the 
orthogonality measure. This is seen by comparing the two integrals at 
the end of the section. 



Consider the matrix 





1 





\ 


Q 





p 







<1 


p 




\ 






■J 



with < p < 1 and q = 1 — p. We look for polynomials Qj (x) such 
that 

Q_i(ar) = 0, Q (x) = 1 
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and if Q(x) denotes the vector 



Q(x) = 



/Qo(x)\ 
Qi(x 

V i / 



we ask that we should have 



FQ(x) = xQ(x). 



The matrix P can be conjugated into a symmetric one and in this 
fashion one can find the explicit expression for these polynomials. 



We have 
Qj(x) = 



i/2 



where Tj and Uj are the Chebyshev polynomials of the first and second 
kind. 

If p > 1/2 we have 



P 



1 — p 



Qn{X)Qm{X)— dX = d r , 



1-X 2 



2(1 -p)n, n = 
2p(l — p)n, n > 1 



while if p < 1/2 we get a new phenomenon, namely the presence of 
point masses in the spectral measure 



p y [" r^™ t 



1-x 2 



dx 



+ (2 - Ap)7T[Q n (l)Q m (l) + Q n {-±)Qm{-l)\ 



2(1 -PK 







2p(l — p)n, n > 1 
From a probabilistic point of view these results are very natural. 

5. The Hahn polynomials, Laplace and Bernoulli 

As has been pointed out before, a limitation of this method is given 
by the sad fact that given the matrix P very seldom can one write down 
the corresponding polynomials and their orthogonality measure. In 
general there is no reason why "physically interesting" Markov chains 
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will give rise to situations where these mathematical objects can be 
found explicitly. 



The example below shows that one can get lucky: there is a very 
old model of the exchange of heat between two bodies going back to 
Laplace and Bernoulli, see [Fj, page 378. It turns out that in this case 
the corresponding orthogonal polynomials can be determined explicitly. 

The Bernoulli-Laplace model for the exchange of heat between two 
bodies consists of two urns, labeled 1 and 2. Initially there are W white 
balls in urn 1 and B black balls in urn 2. The transition mechanism 
is as follows: a ball is picked from each urn and these two balls are 
switched. It is natural to expect that eventually both urns will have a 
nice mixture of white and black balls. 

The state of the system at any time is described by w, defined to 
be the number of white balls in urn 1. It is clear that we have, for 



w = 0,l,...,W 



u^w.w— 1 



W — wW — W 

W B~ 
w B — W + w 

W B 



m wW-w W-wB-W + w 
p = i ! 

W B W B 

Notice that 

Introduce now the dual Hahn polynomials by means of 



■n, —x, x — W — B — 1 
-W,-W 

n = 0, . . . , W: x = 0, . . . , W. 



R n (X(x)) = 3 F 2 



These polynomials depend in general on one more parameter. 

Notice that these are polynomials of degree n in 
X(x) =x(x- W-B-l). 

One has 

/ x(B + W -x+l)\ 
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This means that for each value of x — 0, . . . , W the vector 



\R w {X(x))J 

is an eigenvector of the matrix P with eigenvalue 1 — x( - B+ ^^ x+1 ^ . The 
relevant orthogonality relation is given by 

w 

7ij i?j(A(x))i? :) (A(x))/i(x) 



5ij 



x=0 



with 



w\(-w) x (-w) x (2x-W-B-l) 
(-l)*+i x \(-B) x (x -W-B- l) w+l 



-w 



7T; 



"B) w —j 



j- [w-jy. 

The Karlin-McGregor representation gives 

w 



(f 



7T 



Ri{Kx))RA X (x))e n {x)fi{x) 



x=0 



with e{x) = 1 - 

These results can be used, once again, to get some quantitative 
results on this process. 

Interestingly enough, these polynomials were considered in great 
detail by S. Karlin and J. McGregor, |KMcG2| and used by these au- 
thors in the context of a model in genetics describing fluctuations of 
gene frequency under the influence of mutation and selection. The 
reader will find useful remarks in IDSI. 



6. The classical orthogonal polynomials and the bispectral problem 

The examples discussed above illustrate the following point: quite 
often the orthogonal polynomials that are associated with important 
Markov chains belong to the small class of polynomials usually referred 
to as classical. By this one means that they satisfy not only three term 
recursion relations but that they are also the common eigenfunctions 
of some fixed (usually second order) differential operator. The search 
for polynomials of this kind goes back at least to [B]. In fact this issue 
is even older, and was considered in a paper in 1884 by Routh. See 
[U [R] for a more complete discussion. 
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In the context where both variables are continuous, this problem 
has been raised in |DG| . For a view of some related subjects see [HKJ. 



The reader will find useful material in |AWl IAAR[ Hj. 

7. Matrix valued orthogonal polynomials 

Here we recall a notion due to M. G. Krein, see |Kll IK2| . Given 
a self adjoint positive definite matrix valued smooth weight function 
W(x) with finite moments, we can consider the skew symmetric bilinear 
form defined for any pair of matrix valued polynomial functions P(x) 
and Q(x) by the numerical matrix 

(P,Q) = (P,Q) W = [ P(x)W(x)Q*(x)dx, 

where Q*(x) denotes the conjugate transpose of Q(x). By the usual 
construction this leads to the existence of a sequence of matrix valued 
orthogonal polynomials with non singular leading coefficient. 

Given an orthogonal sequence {P n (x)}n>o one gets by the usual 
argument a three term recursion relation 

(1) xP n (x) = A n P n ^(x) + B n P n (x) + C n P n+1 (x), 

where A n , B n and C n are matrices and the last one is nonsingular. 

We now turn our attention to an important class of orthogonal 
polynomials which we will call classical matrix valued orthogonal poly- 
nomials. Very much as in (Dj IGPT21 IGPT3] we say that the weight 
function is classical if there exists a second order ordinary differential 
operator D with matrix valued polynomial coefficients Aj(x) of degree 
less or equal to j of the form 

d 2 d 

(2) D = A 2 {x)— + A l {x)— + A {x), 

such that for an orthogonal sequence {P n },we have 

(3) DP: = P n *A n , 

where A n is a real valued matrix. This form of the eigenvalue equation 
appears naturally in |GPT1| and differs only superficially with the 



form used in [D], where one uses right handed differential operators. 

During the last few years there has been quite a bit of activity cen- 
tered around an effort to produce families of matrix valued orthogonal 
polynomials that would satisfy differential equations as those above. 
One of the examples that have resulted from this search, see [G], will 
be particularly useful later on. 
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8. Pentadiagonal matrices and matrix valued orthogonal polynomials 

Given a pentadiagonal scalar matrix it is often useful to think of 
it either in its original unblocked form or as being made, let us say, 
of 2 x 2 blocks. These two ways of seeing a matrix, and the fact that 
matrix operations like multiplication can by performed "by blocks" has 
proved very important in the development of fast algorithms. 



In the case of a birth and death process it is useful to think of a 
graph like 




Suppose that we are dealing with a more complicated Markov chain 
in the same probability space, where the elementary transitions can go 
beyond "nearest neighbours". In such a case the graph may look as 
follows 




The matrix P going with the graph above is now pentadiagonal. By 
thinking of it in the manner mentioned above we get a block tridiagonal 
matrix. As an extra bonus, its off-diagonal blocks are triangular. 
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The graph below 




clearly corresponds to a general block tridiagonal matrix, with blocks 
of size 2x2. 

If Pi j denotes the i, j block of P we can generate a sequence of 2 x 2 
matrix valued polynomials Qi(t) by imposing the three term recursion 
of section 8. By using the notation of section 2, we would have 

PQ(x) = xQ(x) 

where the entries of the column vector Q(x) are now 2x2 matrices. 

Proceeding as in the scalar case, this relation can be iterated to 
give 

P n Q(x) = x n Q(x) 

and if we assume the existence of a weight matrix W(x) as in section 
7, with the property 

{Qj,Qj)6ij= [ Qi(x)W(x)Q*(x)dx, 

it is then clear that one can get an expression for the entry of the 
block matrix P n that would look exactly as in the scalar case, namely 

(T) ij (Q j ,Q j ) = J x n Qi(x)W(x)Q*(x)dx. 

Just as in the scalar case, this expression becomes useful when 
we can get our hands on the matrix valued polynomials Qi(x) and 
the orthogonality measure W(x). Notice that we have not discussed 
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conditions on the matrix P to give rise to such a measure. For this issue 
the reader can consult [DPI ID1| and the references in these papers. 

The spectral theory of a scalar double-infinite tridiagonal matrix 
leads naturally to a 2 x 2 semi-infinite matrix. This has been looked 
at in terms of random walks in [P]. In [ILMVJ this work is elaborated 
further to get a formula that could be massaged to look like the right 
hand side of the one above. 



9. An explicit example 

Consider the matrix valued polynomials given by the 3-term recur- 
sion relation 



A n $„_i(x) +B n $ n (x) + C n $ n+1 (x) =t$ n (x),n > 



$_i(a;) = 0, $ (t) = / 



where the entries in A n , B n , C n are given below. 



RANDOM WALKS AND ORTHOGONAL POLYNOMIALS 



15 



Notice that the matrices A n and C n are upper and lower triangular 
respectively. 



in 



A 



12 



A 



B 



21 

n 



n(a + n)(p + 2a + 2n + 3) 



(P + a + 2n + l){(5 + a + 2n + 2)(/3 + 2a + 2n + 1) 
2n((3 + 1) 



(P + 2n + l)(/3 + a + 2n + 2)((3 + 2a + 2n + 1) 
Al 1 := 

22 n(a + n + l)(p + 2n + 3) 



(p + 2n + 1)(P + a + 2n + 2)(p + a + 2n + 3) 
s n . n()9 + n + l)(/? + 2n-l) 
n ' (/5 + 2n + l)(/3 + a + 2n + l) 

(n + l)(/3 + n + 2)(/3 + 2n+l) 
(/3 + 2n + 3)(p + a + 2n + 3) 

2Q3+1) 2 

(/3 + 2n + + 2n + 3)(/3 + 2a + 2n + 3) 

12 2(p + l)(a + p + n + 2) 

(P + 2n + 3){p + a + 2n + 2){p + 2a + 2n + 3) 

2(a + n+ l)(p + 1) 

(P + 2n + 1)(P + a + 2n + 3)(/3 + 2a + 2n + 3) 
n(/3 + n + l)(/3 + 2n + 3) 



S ■= 1 + 

n ' (/3 + 2n + l)(/3 + a + 2n + 2) 

(n+ l)(P + n + 2)(p + 2n + 5) 



+ 



(P + 2n + 3)(p + a + 2n + 4) 
2(/3+ l) 2 



(/3 + 2n + l)(p + 2n + 3)(p + 2a + 2n + 3) 
and finally for the entries of C n we have 

(p + n + 2)(p + 2n + 1)(P + a + n + 2) 



ii 



(/3 + 2n + 3)(P + a + 2n + 2)(/3 + a + 2n + 3) 
2(P+l)(P + n + 2) 



Cn 



a 



22 



(p + 2n + 3)(P + a + 2n + 3)(p + 2a + 2n + 5) 
(p + n + 2)(P + a + n + 3)(p + 2a + 2n + 3) 
(P + a + 2n + 3)(p + a + 2n + A)(p + 2a + 2n + 5) ' 
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If the matrix ^o(x) is given by 

T 1 



= I (f3+2a+3)x _ 2(q+l) 



P+l P+l 

one can see that the polynomials $ n (x) satisfy the orthogonality rela- 



tion 

-l 



§i{x)W{x)$*(x)dx = if i ^ j 



where 



The polynomials $ n (a;) are "classical" in the sense that they are eigen- 
functions of a fixed second order differential operator. More precisely, 
we have 

where 

-n 2 - (a + P + 2)n + a + l + S±L 

-n 2 - (a + (3 + 3)n, 

and T is given by 
T = x(l — x) I — 



+ 



An 



(a+l)(/3+2a+5) _ / p qx 2a+2 , 

/3+2a+3 ^U! -T ^ -T OJX 2a+/3+3 T x 

g±l (a+2)/3+2a 2 +5a+4 / ^ , x 

/3+2a+3 /3+2Q+3 ^T^Tt;A ( 



+ 



g+1 
2 





a + l + 



As mentioned earlier this is the reason why this example has sur- 
faced recently, see [G]. An explicit expression for the polynomials then- 
selves is given in |T2| . corollary 3. 

Now we make two observations 

1) the entries of the corresponding pentadigonal matrix are all 
non-negative. 

2) the sum of the entries on any given row are all equal to 1. 



This allows for an immediate probabilistic interpretation of the pen- 
tadiagonal matrix as the one step transition probability matrix for a 
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Markov chain whose state space could be visualized by the graph given 
below. 




I find it rather remarkable that this example which was produced for 
an entirely different purpose should have this extra property. Finding 
an appropriate combinatorial mechanism, maybe in terms of urns, that 
goes along with this example remains an interesting challenge. 

Two final observations dealing with these state spaces that can be 
analyzed using matrix valued orthogonal polynomials. If we were using 
matrix valued polynomials of size N we would have as state space a 
semiinfinite network consisting of N (instead of two) parallel collection 
of non-negative integers with connections going from each of the N 
states on each vertical rung to every one in the same rung and the 
two neighbouring ones. The examples in |GPT1| give instances of this 
situation with a rather local connection pattern. 

In the case of iV = 2 one could be tempted to paraphrase a well 
known paper and say that "it has not escaped our notice that" some 
of these models could be used to study transport phenomena along a 
DNA segment. 

10. Another example 

Here we consider a different example of matrix valued orthogonal 
polynomials whose block tridiagonal matrix can be seen as a scalar 
pentadiagonal matrix with non-negative elements. In this case the 
sum of the elements in the rows of this scalar matrix is not identically 
one, but this poses no problem in terms of a Karlin-McGregor type 
representation formula for the entries of the powers P n . 
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This example has the important property that the orthogonality 
weight matrix W(x), as well as the polynomials themselves are explic- 
itly known. This is again, a classical situation, see |CG| . 

Consider the block tridiagonal matrix 

/Bo I \ 

A 2 B 2 I 



with 2x2 blocks given as follows 



-I, B n = 0ifn>l 
I if n > 1. 



In this case one can compute explicitly the matrix valued polynomials 
P n given by 

A n P n _x(x) + B n P n (x) + P n+1 (x) = xP n (x), P-i(x) = 0, P (x) = I. 
One gets 

P — ( Un ^ - U n-l(x) 

where U n (x) are the Chebyshev polynomials of the second kind. 



The orthogonality measure is read off from the identity 

i £ Pt(*)-j=^ (I i) = V- 

We get, for n = 0, 1, 2, . . . 

where, as above, 

stands for the i,j block of the matrix P n . 

In this way one can compute the entries of the powers P n with F 
thought of pentadiagonal matrix, namely 
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11. A few more challenges 

We have already pointed out a few challenges raised by our at- 
tempt to extend the Karlin-McGregor representation beyond its orig- 
inal setup. Here we list a few more open problems. The reader will 
undoubtedly come up with many more. 

Is there a natural version of the models introduced by Bernoulli- 
Laplace and by P. and T. Ehrenfest whose solution features matrix 
valued polynomials? 
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Is it possible to modify the simplest Chebyshev type examples in 
|D1| to accomodate cases where some of the blocks in the tridiagonal 
matrix give either absorption or reflection boundary conditions? 

One could consider the emerging class of polynomials of several 
variables and find here interesting instances where the state space is 
higher dimensional. For a systematic study of polynomials in several 
variables one should consult |DX| as well as the work on Macdonald 
polynomials of various kinds, see [M]- A look at the pioneering work 
of Tom Koornwinder, see for instance |Ko| . is always a very good idea. 



After this paper was finished I came up with two independent 
sources of multivariable polynomials of the type alluded to in the pre- 
vious paragraph. One is the series of papers by Hoare and Rahman, 
see [HRT1 IHR21 lHR3l IHR4] . The other one deals with the papers [HI] 
and HXl 



In queueing theory one finds the notion of Quasi-Birth-and-Death 
processes, see |LRll IN] . Within those that are nonhomogeneous one 
could find examples where the general approach advocated here might 
be useful. 
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